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Abstract. A method is presented to search for a hypertorus symmetry axis by the 
alignment of distant objects. This offers greater sensitivity than previously proposed 
object-based methods that rely on accurate true distances. When applied to the catalog 
of objects with z > 1, we find no evidence for a compact dimension. We deduce a lower 
limit to the compact dimension size D > 0.9 of the distance to the cosmic horizon. 
This is consistent with independent constraints from the recent analysis of the WMAP 
microwave background data. 
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1. Introduction 

Large compactified spatial dimensions often feature in modern theories such as String 
Theory, for example in scenarios where the dimensions are all initially small then some 
become large nia. Whether the size of the compact space should be small enough for us 
to detect, however, is much less clear. Observations of the cosmic microwave background 
(CMB) and recent object surveys, have spurred interest in detecting evidence of a 
compact topology by looking for correlations in the positions of distant objects puini, 
or patterns in the CMB lam. The strongest result appears to come from a recent 
analysis of the CMB |H] . From that work it has been deduced that for the hypertorus and 
many other topologies any compact dimension must be larger than 2 cos(25°)r7o ~ 1.8 r^o, 
where rjQ is the distance to the horizon from Earth (half the diameter of the horizon). 
There is, however, uncertainty in the validity of this method due to the thickness of 
the surface of last scattering from which the CMB originates. Some indications were 
previously found of spherical-dodecahedral topology, but these were ultimately not 
found to be statistically signihcant p cni cn eg. The constraints from object based 
methods have been weaker pirns!- Nevertheless, the motivation for a search for a 
hypertorus cosmology has been boosted by possible evidence of a suppression P along 
a common axis in the quadrupole and octopole moments of the CMB. This suppression 
appears signihcant despite the increased cosmic variance at large scales m. While other 
explanations have been put forward |T3l [THIIT7|. there remains the possibility that this 
suppression is due to a sub-horizon scale compact dimension, while taking into account 
the restrictions on dimension ratios introduced by ^2]- In particular we consider a 
compact dimension in hat space as produced by toroidal compactiheation, or as an 
approximation to slightly curved space Hg. For this case, we show that it is possible 
to construct a specialized object-based test with much greater sensitivity than previous 
object-based tests. This provides an independent approach to the CMB methods. 

2. Axisymmetric test 

The test is specihe to cylindrical compactiheation. As discussed below, it can be 
modihed for the other hat-space compactiheations that arise from variations of the 
hypertorus, and even for curved space, but then the upper limit for the detected compact 
dimension is too small to be useful. Figure^shows a plane section through the comoving 
covering space, with the plane containing the compact dimension direction, represented 
by the dotted line. This is a way to visualize the geometry of the compact space, 
formed by imposing cyclic boundary conditions on a non-compact space, euclidean space 
in this case. The dashed lines are on ’copies’ of a surface normal to the direction of 
compactiheation. Arrows represent diherent paths of light rays directed from a distant 
object to the Earth. The copies of an object all make the same angle 0 around the 
compact dimension axis, shown as a dotted line. So, in a perspective projection along 
the axis from Earth, with front and back views superimposed, pairs of images (each 
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Figure 1. Section through the comoving covering space showing light rays from the 
object to Earth 



Figure 2. Perspective view along the compact axis from Earth 


corresponding to two copies of an object in the covering space) are each aligned with 
the center. Figure El illustrates this, with object numbering matching Figure ^ and 
the orientation of 0 indicated. Other object pairs are shown in different colors. This 
observation suggests that we look for an axis direction with an unusually high number 
of such alignments. The same idea appears to carry over to general topologies in 
constantly curved space, from observations using the curvedSpaces softwarej:. However, 
the minimum separation of images that align in this way becomes a multiple of the 
shortest compact dimension, making it weaker than existing limits. 

We shall now describe the test developed here in more detail. The objects used have 
z > 1.0, and are mostly quasars taken from the NASA/IPAC Extragalactic Database § 
(NED). Nearer objects do not help in improving existing constraints on the hypertorus 
dimension. Figure El shows the range and distribution of redshift. A redshift of 5, where 
the data thins out, corresponds to a distance from earth of 0.5 using the currently 

f http: / /WWW. geometrygames.org/ CurvedSpaces/ 

§ http://nedwww.ipac.caltech.edu/ 
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Figure 3. Distribution of redshift among NED objects considered. 


estimated ACDM cosmology. So we might to test a compact dimension np to a distance 
{2){0.5)r]o = r]o. 

The Earth’s pecnliar velocity, its velocity relative to comoving space and the CMB, 
(3 ~ 0.00123, canses a maximnm angnlar error 0.1° and maximnm redshift error 
~ 0.01 [Tin 1^ . This aberration is hrst removed from the snrvey data by applying a 
boost —j3 to incoming light rays. An aberrated line of sight vector h' is deaberrated EH 
to a line of sight vector 


n' + [(7 - l)(3.n' - 7 / 3]/3 
7(1 - /3.n') 


7 = (1-/3^)-^/^ 


( 1 ) 


Where (3 denotes (3 normalized. Similarly, the observed redshift of an object z' 
deaberrates to 


+^1) 

7(1 - (3.n') 


- 1 


( 2 ) 


The error remaining dne to the motion of the Earth relative to the Snn, 30 km/s, and 
the motion at the Earth’s snrface dne to rotation, 0.5 km/s, is from |Tn], (3 radians 
~ 0.006°. In order to detect a compact dimension of size np to the horizon, we look 
for image matches on opposite sides of the sky for each possible axis of the compact 
dimension. The alignment check is made more efficient by maintaining a list of objects 
ordered by angle, 0, abont the search axis. Alignments are then fonnd by rnnning 
throngh the list only once. After shifting the axis, the list is reordered slightly nsing 
a fast sorting algorithm. Dne to the hnite angnlar resolntion of the search Aa, each 
object is given an allowed range of 0 valnes [0 — A(j), 0 + Acp] where A0 = Aa/ sin^ 
and 6 is the angle of the object away from the axis, as shown in FignreEl When the 0 
ranges of two objects on opposite sides of the axis overlap, a candidate match has been 
fonnd. 
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Figure 4. 



Time, Gyr 



z, redshif t 


Figure 5. Look-back time functions for = 0.17,0.3, 0.42(upper). 


To reduce background noise from chance alignments, further tests are performed 
once an alignment has been found. Quasars are not expected to live longer than 10® 
years ji], so images that differ in look-back time by more than this are discarded. The 
look-back time is determined from the redshift. We assume a flat ACDM model with 
0.17 < Dm < 0.42. Variations in the look-back time caused by uncertainty in the 
cosmological model can be tolerated, since this has a minor effect on the rejection of 
pairs due to look-back time difference. Hence, a flat Dm = 0.3 ACDM model is used. 
The similarity of the look-back time differences can be seen in Figure El which shows 
the curves for the extremes of Dm enclosing the Dm = 0.3 curve. A function is used to 
accurately approximate the exact integral expression for look-back in a flat ACDM |22] 
universe given by 

t{z,Dm,D/^) = ^ [ dz{l +z)~^{Dm{l +zf+ Dk)~^/'^ . (3) 

Hq jo 

Taking the age of the universe as 13.7 Gyr [221, approximate look-back formula used 
with Dm = 0.3 is 

t{z) = 23.8(^(1 - (1 + z)-^/^)) - 2.3 Gyr . 


(4) 
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Figure 6. Look-back time function approximation error. 


This look-back time yields look-back time differences accurate to within 0.1 Gyr 
up to ^ = 8 for the given range of cosmologies, as shown in Figure IHl The look-back 
hlter also has the very positive effect of greatly reducing the angular error caused by 
the peculiar velocity of the object to a value which is well below other errors. This 
is because matching images at the same distance will originate at the same time, and 
hence from the same place in the compact space, regardless of the peculiar velocity. 
Finally, we notice that the perpendicular (covariant) distance from each image to the 
axis is the same, giving the constraint sin(6'i)d(^i) = sin{ 62 )d{z 2 )■ Allowing for angular 
and redshift errors, this translates into requiring an overlap of the two intervals, 

[sin( 6 'i — Aa)d{zi — A 2 ;),sin( 6 'i -|- Aa)d{zi A^)], 

[sin( 6*2 — Aa)d{z 2 — Az), sin( 6*2 -|- Aa)d{z 2 + Az)]. 

We choose Aa = 0.01° consistent with the residual errors from aberration and 
peculiar velocity, and A 2 ; = 0.01 consistent with the observational error of many 
spectrographic redshift values. The test is sensitive to the cosmology determining 
the covariant distance function d{z). Therefore it is repeated for a range of ACDM 
models about flm = 0.3. d{z) is calculated using an interpolated lookup table for each 
cosmology. The tables are calculated using the following formula for distance 

d{z) = [ dz{{lzf{l QrnZ) - flAz{2z))~^^‘^) , (5) 

JO 

normalized to 1 at the horizon. 

We also found that one hnal hlter is necessary to exclude false peaks which result 
from tight clusters of correlated objects. These exist due to duplication of objects in the 
database and the high concentrations of objects found in pencil-beam surveys. These 
are excluded by organizing pairs into groups that are separated by a minimum angle of 
0.5°. The angle was chosen because it is sufficient to remove the clustering effect, while 
giving a very small chance of reducing the numbers of genuine matching pairs. 

Once a pair of images has passed all tests, the compact dimension size is estimated 
as D = di cos 6^1 -|- ^2 cos 02 , (cf. Figure Ej). The alignment is used to increment the 
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bin corresponding to D in a histogram corresponding to 0^- The image references 
are also added to the bin, for possible later query. Once the search over (j) is 
complete the histograms are analyzed to see if any bins contain more alignments than 
the existing record values. If so, the record values are updated, including image 
references. The search then proceeds to the next candidate axis. The HEALPix || 
function pix2ang_ring() |21] is used to conveniently generate a uniformly spaced set 
of directions covering half the sky, and therefore all axes. The ring ordering ensures that 
all the angle increments are small, minimizing reorderings. The HEALPix resolution is 
chosen to be = 5000 giving a pixel spacing Aa ~ 0.01° consistent with the errors 
from the aberration and peculiar velocity. 

3. Simulated alignments 

As a test of the method and to calibrate the bin sizes we added simulated aligned 
objects to the database. The D bin and flm spacing are made just large enough so that 
a simulated set of aligned images, of look-back difference less than 1 Gyr, with maximum 
expected noise added and maximum misalignment between the values, can be fully 
detected in one bin. The 26 cosmologies were tested were with 0.17 < Qm < 0.42 in 
steps of 0.01. The distance D, in units of distance to the horizon, was binned from 0.2 to 
1.2 in steps of 0.002. The bins were overlapping to guarantee maximum detection rate. 
Figure [3 shows the output of record alignments with 20 simulated object pairs added 
to the objects data, and D = 1.0, for the closest cosmology tested. The additional lines 
are signihcance contours described in the next section. 

4. Results 

The test was applied to a catalog of 40000 objects with 2; > 1.0 obtained from the NED 
database. The search space was divided evenly between 20 1.4 GHz Linux workstations, 
and ran for 5 days. There was no clear evidence of a peak. To illustrate this. Figure 
IHl shows the maximum number of alignments found in each D bin for 12^ = 0.3, once 
all the searches are combined. The contours are for a, 2a, Sa signihcance as described 
below. The cost of searching prohibits a detailed characterization of the noise statistics 
by repeated searches over decorrelated data sets. We expect only one anomalous peak, 
or possibly a few if we consider a hypertorus with more than one dimension inside 
the horizon, but all the plots show very similar local noise distributions. To quantify 
this more precisely we need to estimate the noise statistics from the result set and 
calculate how big a peak must be to be signihcant. We do this by looking at the 
huctuation in alignments in each local region of D. The task is complicated by the 
number of alignments being a discrete variable. However, looking over the whole range 
of cosmologies we hnd that the local standard deviation in the noise, ai, varies from 
~ 0.7 in the central regions to ~ 0.3 at the sides. The distribution for the maximum 

II http://www.eso.org/science/healpix/ 
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Figure 7. Test of 20 simulated pairs with D = 1.0, plus observed data. 


of a set of Poisson variables, which approximates the record alignments, falls off faster 
than a normal distribution. Therefore, we can make a conservative estimate of the 
error contour above which peaks are signihcant by using a normal model for the local 
noise distribution with Let the contour be at na^ above the local average of the 
noise. Then the probability of hnding any peaks above the error contours for any of the 
cosmologies is 1 — {1 — ErruT^n))^ ^ NErrRxin), where N =(num. cosmologies)(num. 
D-bins)= (26)(500) = 13000 is the total number of bins in all cosmologies, and Erruxif) 
is the right tail normal error function. To make this probability correspond to an 
overall 3a signihcance we set it to 0.003, obtaining n ~ 5.0 and nai ~ 3.5 in the 
central region. The contours corresponding to cr, 2a, 3ct signihcance are shown in Figure 
0 To summarize, the chance of hnding any peaks above the 3a contour, for any of 
the cosmologies, is 0.003. We have given conditions for a statistically signihcant event. 
Next we need to show that the conditions will be met if a compact dimension actually 
exists. 

5. Sensitivity 

The question remains as to whether the test just described is sensitive enough to detect 
a compact dimension, and can we identify parameters in the catalog used that ahect 
sensitivity. To answer this we hrst determined how many alignments would be generated 
along a given axis, for the full range of compact dimension sizes D, given observations 
on one half of the axis from the NED data, and assuming every observed image has a 
matching observable partner if they can satisfy the look-back time diherence criteria. 
We dehne completeness as the fraction of those observed images that could have a 
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D, compact dimension size 


Figure 8. Results of search compact dimension search for = 0.3 


matching partner satisfying look-back difference criteria, which actually do have an 
observed partner. Then the last assumption is equivalent to saying the completeness 
is 1, or we also say the data is complete. The result for a typical axis with flm = 0.3 
is shown in Figure IHl Figure fTTII shows the ratio of the complete alignments to the 3a 
error contour for fl^n = 0.3, which is the degree of incompleteness, the reciprocal of 
completeness, that can be tolerated while still making a positive detection. In other 
words from the initial calculation assuming completeness, we are calculating the lowest 
levels of completeness at which a detection can still be made for each D. For much of the 
D range this is much higher than the incompleteness tolerated in [Hj. In their work it is 
claimed that a toroidal topology can be detected with rejection of 90% or completeness 
of ~ 0.1 in the current catalog. However, this rejection rate includes rejections due to 
look-back time difference, whereas in our treatment the completeness quoted is for the 
original data set. In practice look-back difference rejection would reduce the overall 
completeness well below their rejection threshold, so that even with a complete original 
set their test would not be sensitive enough. 

The next question is can we estimate the incompleteness of our data? One way 
to approach this is to look at the density of objects in detailed pencil surveys as an 
indication of total number of objects we can detect, and compare this with the average 
across the sky. Doing this for a few pencil surveys gives a density ratio of ~ 10^ implying 
that the data is too incomplete to detect a compact dimension. However this probably 
overestimates the incompleteness because the pencil surveys include objects of much 
lesser magnitude than the remaining surveys, and we expect matched objects to be of 
similar absolute magnitude, because of similar look-back times. So, it is possible that 
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Figure 9. Maximum alignments along one axis for the NED catalog made complete. 



Figure 10. Maximum incompleteness for which a significant detection can be made. 
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the survey is complete enough to make detection possible over a wide range of D: From 
Figure Uni a completeness of 0.02 would rule out a compact dimension up to 0.9 r/o- 

The dependence of sensitivity on completeness can be modeled as follows. Suppose 
a complete observable set contains n matched pairs, and the fraction of images actually 
observed or completeness is p. Then np^ matched pairs will be observed on average. The 
number of pairs observed by chance can be checked empirically to be of the form z/p^, 
for a constant v. This can be explained as the maximum of a set poisson distributions, 
which is proportional to the mean of a single poisson distribution describing the number 
of alignments in a given direction. The mean is clearly proportional to p^ since an 
alignment occurs when events on two linear poisson processes intersect. The upshot 
is, perhaps surprisingly, that the ratio of signal to noise njv is independent of p, the 
completeness, except that additionally we require p sufficient to ensure that the signal 
is separated from the fluctuations in the noise. This means that we are immune to 
incompleteness. The calculation of incompleteness previously was used to establish the 
range of D to which we are sensitive, and this range will remain valid for higher levels 
of completeness. It is therefore unnecessary to process more than a certain number of 
images, a random selection suffices. The only way to improve sensitivity is by choosing 
more selection criteria that can raise njv. This will reduce the detected signal and noise, 
but may require more objects so that they can be separated. 

6. Summary 

The approach described is more sensitive than previous object-based tests applied to 
the hypertorus, because it does not rely on accurate true distance determination. Al¬ 
though there is uncertainty over the degree of completeness among the observations, the 
results presented here suggest that a hypertorus cosmology is ruled out for D < 0.9 r^o- 
Selection could be improved by matching magnitude and spectral profiles, using greater 
an enlarged catalog. Deeper observations going beyond the current 0.5 limit will 
extend the range of D tested. 
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